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Abstract 

We compute finite volume effects relevant for lattice QCD simulations using background fields. 
Focusing on constant electric and magnetic fields on a periodic lattice, we determine volume cor- 
rections to pion two-point functions using chiral perturbation theory. Among such corrections are 
the finite volume shifts to the electric and magnetic polarizabilities, which are numerically shown 
to be non-negligible. We additionally find that all terms in the single-particle effective action can 
be renormalized by infrared effects. This includes Born couplings to the pion current and total 
charge-squared, which can be renormalized due to the nature of gauge invariance on a torus. 
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I. OVERVIEW 



Lattice gauge theory simulations are making dramatic progress towards addressing quan- 
titatively the non-perturbative dynamics of quarks and gluons Despite its successes, 
lattice QCD suffers from systematic errors due to a number of sources. Understanding and 
accounting for these systematic errors is necessary for accurate determination of hadronic 
observables. In particular, as the lattice pion masses approach the physical pion mass on 
a fixed-size lattice, the effect of finite volume will become pronounced. This is increasingly 
important to note as lattice calculations are in reach of the physical point [2|. Certain ob- 
servables are more susceptible to volume corrections; thus, a practical application of field 
theories in finite volume is the study of finite-size scaling of observables calculated using 
lattice QCD simulations. 

The long-range physics of low-energy QCD is a consequence of the lightest degrees of 
freedom, the pseudo-scalar pions. In finite volume, hadronic observables are altered due to 
modified pion dynamics arising from boundary conditions imposed on the underlying lat- 
tice action. These finite-size effects can be treated systematically using chiral perturbation 
theory [1, 0, [5|. In this theory, the light pions emerge as a consequence of spontaneous 
chiral symmetry breaking, and the smallness of the up and down quark masses compared 
to the QCD scale. For QCD in the presence of external sources, low-energy hadronic prop- 
erties are encoded in multipole moments, radii, and polarizabilities. These quantities are 
accessible with lattice QCD correlation functions using two different techniques: the source 
insertion method and the classical external field method. The source insertion method does 
not directly probe such observables: lattice data additionally require momentum extrapo- 
lation. Furthermore, at finite volume the multipole decomposition into momentum-transfer 
dependent form factors is no longer valid. Finite volume studies of hadronic current matrix 
elements at non-vanishing momentum transfer and their generalizations have been studied 
using chiral perturbation theory jl, 0, 0, 0] • 

In this work, we address finite volume effects relevant for QCD simulations in background 
fields, see e.g. [lol, IHI, [H, M, H H H H H H H [2l|. Focusing on constant electric 



and magnetic fields, we use chiral perturbation theory to calculate pion two-point functions 
in the presence of classical external fields. These results can be utilized to analyze lattice 
correlation functions calculated in background field simulations. Our computations demon- 
strate that the nature of gauge invariance on a torus leads to finite volume artifacts. A 
striking manifestation of this is the infrared renormalization of pion currents and charge- 
squared couplings. Identical effects have been found from momentum space calculations of 



current matrix elements |22l . |23| . Consequently two-point functions for both charged and 
neutral particles are modified at finite volume; for example, neutral particles propagating in 
external fields no longer have correlation functions with a simple exponential falloff at long 
times. On general grounds, one expects the finite volume renormalization of all terms in the 
effective action for a particle in a background field. 

Our presentation is organized as follows. First in Section [Tll we consider the case of 
pions propagating in a background electric field. Here the field is treated in Euclidean 
space. Corrections to the infinite volume single-particle effective action are determined 
using spatially periodic boundary conditions. In Section llllt we consider the case of pions 
propagating in a background magnetic field. The effects of finite volume are determined on 
the single-particle effective action for this case as well. For both electric and magnetic cases, 
we provide expressions for the volume effects under the lattice approximation of omitting 
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disconnected quark contractions. Finally in Section IIVI we summarize our findings, make 
important clarifications, and discuss further applications of our work. 



II. MESONS IN FINITE VOLUME: ELECTRIC CASE 

Pions propagating in a background electric field at finite volume can be handled us- 



ing Schwinger's proper time method [2J]. We choose to implement the field with a time- 
dependent gauge potential 

= (0,0,-^X4,0), (1) 

which corresponds to a Euclidean electric field, in the z-direction. As we are interested 
in quantities perturbative in the field strength, we can perform the analytic continuation to 
Minkowski space trivially j25|. The lattice simulations, moreover, are performed in Euclidean 
space for which the choice in Eq. ([1]) is natural. Throughout, we work with theories defined 
on a spatially periodic torus of infinite time extent. The Euclidean space action has the 
form 

dxi / dxC{x). (2) 

oo Jo 

The length of each spatial direction is given by L. As the field is implemented by a time- 
dependent gauge potential, there is no conflict between Eq. ([1]) and the spatial periodicity 
of the underlying lattice action.^ 

The low-energy dynamics of pions is described by chiral perturbation theory [30], in 
which pions are the Goldstone modes emerging from spontaneous chiral symmetry breaking.^ 
These modes are non-linearly realized in E, where S = exp{2i(j) / f) , and is an SU{2) matrix 
of pion fields. The dimensionful parameter / can be identified with the chiral limit value 
for the pion decay constant, and / = 130 MeV in our normalization. The chiral Lagrangian 
appears as 

= 4 < D^^^D^E >-^< mQ(St + S) >, (3) 
o Z 

where we have used a bracketed notation to denote flavor traces: < A >= tr(yl). The quark 
mass matrix is mq, where rriQ = dia.g{mu, md), and the action of the covariant derivative is 
specified by 

D^^ = d,J: + teA^[Q,j:], (4) 

where Q is the quark electric charge matrix, Q = diag(2/3, — 1/3). In writing this La- 
grangian, we have included only the lowest-order terms. These terms scale as 0{e'^), where 
e is a small dimensionless number. We assume the power counting: k"^ / ~ m^/A^ ~ e^. 

Here A^ is the chiral symmetry breaking scale, A^ = 2-\/2vr/, is a typical momentum. 



^ In actual lattice simulations, the length of the time direction is longer than the spatial directions, but 
obviously not infinite. To arrive at a periodic action in this case, the field strength must be quan tized in 
the form q£ = 2im/ j3L, where (3 is the temporal extent of the lattice, and n is an integer [26l. b?!. [28 . 29 1. 



We assume this quantization condition has been met, but that /3 > 2L so that effects from the finite 
temporal extent of the lattice are comparatively quite small. 
^ At finite volume, we must also specify that the volume is large enough to prevent the pion zero modes from 
conspiring to restore chiral symmetry. With this assumption, we work in the so-called p-regime 0, 0, S] • 
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FIG. 1: Pion propagator in strong external fields. Depicted are the Born couplings to the charged 
pion which have been summed to arrive at Eq. 




FIG. 2: Diagrams for the pion two-point function in a non-perturbative external field. The dashed 
lines represent pion propagation in the background field. The wiggly lines terminating in crosses 
represent couplings of the vertices to the background field. 



and is the pion mass. The leading-order pion mass can be determined by expanding the 
Lagrangian in Eq. ([3]) to quadratic order in the meson fields. Using the strong isospin limit, 
rUu = rrid = m, we have /^m^ = 4Am. 

At leading order, the neutral pion propagator maintains a Klein-Gordon form due to the 
absence of charge couplings to the external field. For the charged pions, there are additional 
interactions. Although we shall be interested in small fields, i.e. eF^u/m^ ~ e^, it is actually 
easiest first to treat this ratio as 0{e^), and then expand quantities to the desired order in 
the external field strength. Summing the Born couplings to all orders, see Figure [1], produces 
the propagator of the positively charged pion 



X 



2L3 



X4 



e8 



.0 



(5) 



where E'j,^ = k\ + + m^, and the periodic momentum modes are given in terms of 
triplets of integers n, by the relation k = 2nn/L. This propagator depends on the quantum 
mechanical harmonic oscillator propagator evolved in proper-time s, 



{X',s\X,0) 



eS 



27r sinh e£s 



exp 



eS 



2 sinh eSs 



[{X'^ + X^) coshers - 2X'X] | . 



(6) 



To compute finite volume corrections to pion two-point functions in background fields, we 
use the charged pion propagator above in loop diagrams that arise in the effective theory, see 
Figure [2J The neutral pion contributions to such loop diagrams only lead to field independent 
modifications. These modifications are part of the finite volume corrections to the pion mass. 
Computation of charged pion loop diagrams with Eq. leads to ultraviolet divergences, 
which we regulate using a simple subtraction 



dk 



(7) 



that produces only the infrared corrections we seek. 

It is instructive to consider the simplest loop contribution explicitly. The four-pion vertex 
arising from the mass term in Eq. ([3]) makes contributions to diagram C depicted in the 
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figure. This mass-insertion tadpole diagram depends on the self-contracted propagator. 
Denoting the space-time coordinate of the vertex as y^, we have the ultraviolet regulated 
pion self-contraction 





,4-^,0 e-^-^x (8) 



fc 

(9) 

where in the second line we have series expanded in powers of the field strength S, keeping 
only the second-order term. The zeroth-order term in S combines with the corrections from 
neutral pion loops to give the finite volume pion mass shift. We assume that the background 
field strengths employed on the lattice are sufficiently small so that the fourth order and 
higher-order terms can be neglected. The sums in Eq. are standard in the study of finite 
volume modifications to meson properties. Using the catalog of finite volume functions 



of [22(], we can write the finite volume modification at 0{S'^) as 



h[Diy,y)] = ^e^£:^ !^yl [J5/2{m^) - h/2im^)] + ^J7/2{m^) - ]j^h/2{m^)^ . (10) 

The last two terms give rise to a field-dependent shift in the effective mass squared. These 
are thus contributions to the finite volume electric polarizability. 

Notice the remaining dependence on the vertex time, y^, in the first terms of Eq. ( ITOl) . In 
infinite volume, such dependence disappears owing to an allowed shift of the integrand: 

"OO POO 

dks e-^^^'-^^y^^' = / dks e-"^3 . 

-OO J — OO 

At finite volume, discrete translational invar iance is not enough to remove the time- 
dependence.^ This time dependence leads to an infrared renormalization of the single-pion 
effective action, as we now demonstrate. Consider the full pion two-point function G{x',x). 
Because we have discrete translational invariance in space, we can choose to locate the pion 
source at a; = 0. Furthermore, because the three momentum is a good quantum number, 
we can project onto a final-state pion at rest. This leads us to 

G{t,0) = [ dxG{x,t;O,0), (11) 
Jo 

and similarly for the perturbative propagator, D{t,0). This notation allows us to focus on 
the time-dependence. A perturbative contribution to G{t, 0) with time-dependence as in 
Eq. ffTOj) has the form 

Git,0) = Dit,0)-A£^ j dy^Dit,y^)ylDiy,,0). (12) 



^ There is one special case: for a discrete torus with field strengths quantized in the form e£ = 27rn/L 
the time-dependence can be shifted away by reindexing the sum on momentum modes. As these field 
strengths are prohibitively large, however, we will not separately address this case for which we would 
need an entirely different treatment. 
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The coefficient A arises from the one-loop graphs; and, as such, can be treated as a pertur- 
bation. Using a quantum mechanical notation, we have 



{t\G\0) = {t\D\0)-A£' j dy,{t\D\y,)yl{y,\D\0) 

= {t\ [D - AS'DT'D] |0) = {t\D^^-^^^\0) + . . . , (13) 
where T\t) = t\t). Applying the operator identity 

^ D-'[D-' + A\-\ (14) 



1 + AD 

we see that the inverse propagator G^^ has been additively renormalized 

G-' = + A£^T^. (15) 

For the case of charged and neutral pions projected onto vanishing three-momentum, the 
inverse perturbative propagator has the form 

D-' = --^ + ml + Q'£^T^ (16) 

where Q is the pion charge. Hence the perturbation by A corresponds to an infrared renor- 
malization of the field-squared coupling: — * + A. This is the coordinate space analog 
of the infrared renormalization of the zero-frequency Compton scattering tensor, originally 
found by j22|. Gauge invariance is maintained, however, because Ward identities no longer 
protect zero- frequency scattering [23|. A gauge invariant operator that gives rise to a shift 
in the charge-squared coupling of the neutral pion, for example, is 

£= -l^^(L)7rVVF(-) (17) 
The W^/""* is defined in terms of Wilson lines, wj; ^ = ^ — , where Wi is given by 

Wi = exp (i [ dxiAix)] . (18) 



The Einstein summation convention has been suspended in this equation. Because the space 
is compact, the Wilson line Wi is actually a gauge invariant loop. The coefficient of this 
operator, A{L), runs to zero exponentially fast as the volume is taken to infinity. Expanding 
the operator out to second order in the background field, we have 

£ = l^(L)7rVfV + 0(f^), (19) 

which has precisely the form needed above to describe renormalization of the charged-squared 
coupling. 

The remaining contributions to the finite volume two-point function arise from expanding 
out the kinetic term in Eq. ([3]) to fourth order in the pion fields. The two gauge covariant 
derivatives generate all three diagrams shown in Figure [2l The gauge part of the derivative 
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leads to an explicit power of the external field at the vertex. Accordingly we can evaluate 
each diagram by using the propagator in Eq. ([5]) in the appropriate contractions, and expand 
the result to second order in the field strength. Unlike the infinite volume calculation, there 
are no particularly useful cancellations to be aware of. Focusing on the sector with vanishing 
three-momentum, the resulting single-particle effective action to 0{£'^) has the form 



/oo 
dt Tx\t) 
■oo 



n{t). (20) 



Here we have used a^; to denote the electric polarizability which gives a positive contribution 
to the effective mass due to our Euclidean space treatment. Each of the renormalized 
couplings is a sum of the infinite volume piece plus a finite volume correction, 

ml{L) = ml + 6L[ml], 
aE{L) = aE + 5l[oie\, 

Q\L) = Q' + dLlQ']. (21) 

With the exception of the pion mass, which we do not consider here, the finite volume 
corrections are isospin dependent. For the neutral pion, we find 

SlIQIo] = -^mt/5/2(m.), (22) 

= ^Z^^' [2/5/2 (m.) + Jr/2(m.)] , (23) 

while for the charged pions, we have 

hlQU = [2/i/2(m.) + mlly.im^)] , (24) 

SL[a]t] = ,^^mlh,2{m^). (25) 

Here a/.^. = e^/47r is the fine-structure constant. For both charged and neutral pions, the 
values of found here agree with those deduced from the zero- frequency Compton 

scattering tensor Focusing on the charged pion electric polarizability, we recall that 

delicate cancelations lead to only neutral pion loop contributions in infinite volume. These 
contributions do not depend on the background field; thus, there are no one-loop effects 
contributing to the polarizability. In finite volume, however, such cancelations are spoiled 
due to the lack of SO (4) invariance. Consequently the value of 6i\a^^ ] is non- vanishing. 

To investigate the size of the finite volume corrections, we consider the relative change in 
each of the quantities. For the charged pion, the relative change in charge squared is just the 
finite volume effect, A(5^± = (5l[(5^±], because in infinite volume = 1. For the neutral 
pion, we must define the relative change, since the infinite volume charge- squared coupling 
is zero. We choose the simplest possible definition, AQ^o = 5l[Q^o]- The charge-squared 
volume effects are shown in Figure [3l To consider the finite volume effects for the electric 
polarizabilities, we must recall the one-loop infinite volume results (31I . 32 . 33l |: 



='"^--r;;^' =~ 3(4./)^m/ 
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FIG. 3: Relative change in charge-squared couphngs. Plotted versus L for a few values of the pion 
mass are the finite volume effects AQ^ for both neutral and charged pions. 
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FIG. 4: Relative change in electric polarizabilities. Plotted versus L for a few values of the pion 
mass are the finite volume effects Aa^; for both neutral and charged pions. 



Because both of these are non-vanishing, we can define the relative change due to volume 
effects in the ordinary manner, Aa^; = ■ To determine the low-energy constants 



required for the charged pion polarizability, we use [34 



ag + aio = ^^/^/ fv^ (27) 



where the form factors are given by = 0.0115 [35j, and fv = 0.0262 ^]. The latter 
value is consistent with the determination of fv assuming conservation of vector current, 
fv = 0.0259. Notice the combination ag + aio is renormalization scale independent. The 
finite volume effect on the electric polarizabilities is show in Figure HI which depicts non- 
negligible corrections. The relative change in neutral pion polarizability is more sensitive to 
volume effects than the charged pion. One reason for this difference is the scaling by the 
infinite volume result. The infinite volume polarizability is a factor of five times smaller for 
the neutral pion compared to the charged pion, while the finite volume modifications are 
roughly the same size. 

The neutral pion polarizability is already challenging to calculate on the lattice due 
to disconnected quark contractions (which have so far been neglected in simulations). In 
infinite volume, the one-loop expression for the neutral pion polarizability stems entirely 
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from such disconnected quark contractions |8| . The electric polarizabihty of the charged pion 
can be determined from the lattice. Reports of a preliminary calculation have been given 



m 



20]. At this stage, contributions from background field couphngs to the sea quarks have 
been omitted due to computational restrictions. To address the omission of disconnected 
quark contractions for neutral and charged pions, we turn to a partially quenched theory to 



differentiate between valence and sea quarks [37|, |38(]. We perform the partially quenched 
finite volume analysis with degenerate valence and sea quarks using a partially quenched 
charge matrix that distinguishes between valence and sea quark charges |39|, |40| . The latter 
are set to zero. Additionally for the neutral pion, we perform the computation for only the 
connected part of the two-point function. The results are as follows: the charge-squared 
couplings are renormalized by 

*4«y«Er" = o. = i4Qy, (28) 

where the subscript QEM denotes that we have quenched the electromagnetic charges, and 
the superscript connected denotes that only the connected part of the correlator has been 
retained. The shifts in polarizabilities are given by 

Sda-EYQEu''"' = h[al%EM = ^^mlh„{m^), (29) 

and are thus of comparable size to those with unquenched charges. Finite volume results 
for pion electric polarizabilities suggest that as the pion mass is brought down on fixed sized 
lattices, it will become increasingly challenging to isolate the infinite volume physics. A 
thorough study of the pion mass and volume dependence of lattice data will be required to 
extract pion polarizabilities. 



III. MESONS IN FINITE VOLUME: MAGNETIC CASE 

Now we investigate the finite volume modifications to pion two-point functions for the 
case of a magnetic field. To be concrete, we implement the background field with the gauge 
potential 

= (-5x2,0,0,0), (30) 

which corresponds to a magnetic field, i?, in the z-direction. On a spatial torus, the imple- 
mentation of Eq. fl5U]) generally breaks the periodicity of the lattice. To maintain discrete 
translational invariance of the lattice action, the field strength must be quantized in the 
form 0, H H [13 

qB = (31) 

where n is an integer. Provided this condition is met, we can expand hadronic field operators 
in periodic momentum modes, and treat them in a constant background field. 

The Schwinger proper-time method cannot be directly utilized for the magnetic case 
because the ?/-component of momentum is missing in the resummed propagator. Pions 
propagating in a background magnetic field at finite volume, however, can be handled us- 
ing the free particle propagator in coordinate space. We return to Eq. ^ and treat the 
external field as a small perturbation, B/m\ ~ e^. Thus to leading order, the external field 
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FIG. 5: Diagrams for pion two-point functions in a perturbative external field. The straight lines 
represent free pion propagation. The wiggly lines terminating in crosses represent interactions with 
the background field. 

dependence can be neglected, and we have the coordinate space propagator 

1 , , r rlh. pik-{x' -x) ik4{x'^-X4) 

n -in 

for both charged and neutral pions. The spatial momentum modes are again quantized in 
the form k = 2TTn/L. The remaining field-dependent Born terms in Eq. must then 
be treated as coordinate-space interactions in perturbation theory. There are two types 
of interactions: a two-pion, linear background field coupling; and, a two-pion, quadratic 
background field coupling. 

To calculate one-loop corrections to pion two-point functions, we expand the Lagrangian 
in Eq. ([3]) to fourth order in the pion fields. We then generate all possible one-loop diagrams 
up to quadratic order in the external field B. Due to symmetry, diagrams linear in the 
external field vanish.^ The non- vanishing diagrams are shown in Figure [51 It is particularly 
illustrative to consider a simple contribution explicitly. Consider the four-pion vertex with 
mass insertion derived from Eq. Q. Using this vertex along with the quadratic background 
field coupling, we arrive at a contribution to diagram C in the figure. Let us denote the 
location of the four-pion vertex as and the location of the quadratic field insertion as 2;^. 
Then the correction to the propagator, SG{x',x), from this diagram is of the form 

/oo 
dy, / dyD{x',y)Aiy)D{y,x), (33) 
-00 Jo 

where, up to overall constants, 

/oo rL 
dz, / dzD{y,z)zlB^D{z,y). (34) 
-00 Jo 



^ This is true when the action is projected onto vanishing x-component of momentum, k\ — 0. For a 
general momentum, the infinite volume action contains a term proportional to 2iBx2ki, which acquires 



infrared renormahzation in finite volume. This corresponds to screening of charged particle currents 23[. 
Charge conjugation invariance forbids this effect for the neutral pion. 
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Inserting the perturbative propagator, Eq. (132|) . into the expression for A{y), yields three 
trivial integrations that can be performed: those over zi, Z3, and Z4. Clearly the Z2 integra- 
tion is non-trivial given the coordinate dependence introduced by the gauge potential. After 
performing trivial integrations, we have 



A(2/) 



dZn 



E 



^2 ^ik2(y2-Z2) ^ik'^(z2-y2) 



27r + kl + ml] [kl + k'i + kl + kl + ml] 



E 



-i(k2—k'2)Z2 



92 



oik2y2 



dk4 

'^6:^1^1 + 4' + kl + kl + ml] \dkl + kl + ml 



-ik'2y2 



(35) 



with k = 27rn/L, and fcg = 2iim/L. We then differentiate, perform the integrals over Z2 and 
^2 that have been rendered trivial, and drop terms that vanish by parity. Regulating the 
sum over momentum modes, we are left with 



B' 



B' 




vl 



+ 



k\ + m. 



214 



-^1/2 ^3/2 ("^^) 



3 5 

-4/2 ("^tt) + —^7/2 ("^7, 



(36) 



The last two terms are finite volume corrections to the magnetic polarizability, while the first 
term is a finite volume correction to the charge-squared coupling. The argument exposing 
this fact follows, mutatis mutandis, from that for the electric case. 

The remaining one-loop contributions to the pion two-point function can be evaluated 
similarly. The net result is that the pions are described by an effective action of the form 



S, 



eS 



dt 



dy TrHy,t) 



dt^ 



8^ 



ml{L) - Anm^/SMiL) B'' + Q'\L) e'B^ 



7r(l/,t), 
(37) 



up to 0{B^) corrections. In writing this effective action, we have projected onto the sector 
with kl = ks = 0. The modification to the charge-squared term of the action is allowed 
by gauge invariance due to exactly the same operators as in the electric case. The single- 
particle effective theory consequently demands that the charge-squared couplings be identical 
to those determined above [23] . The results of the explicit finite volume calculations indeed 
show that the Q^{L) in Eq. fl37|) are identical to those in Eq. fl20|) for both 7r° and tt^. The 
same is additionally true for the case where disconnected quark contractions are omitted. 

The term in the effective action with coefficient I3m{L) denotes the magnetic polarizability 
contribution. This coefficient is a sum of the infinite volume polarizability, and the finite 
volume modification 

(3m{L) = (3m + hi^M]. (38) 

The former is given hy [3m = ~CiE for charged and neutral pions, due to the helicity structure 
of the one-loop diagrams, while the latter are given by 



h[(3ln 



M 



24/2m^ 
12/%^ 



I [2/5/2 (m^) + 7Kg/2{m^ 



m. 



m^J7/2(m^ 



(39) 
(40) 



11 




2.5 3.0 3.5 4.0 4.5 5.0 5.5 6.0 



0.0 



2.5 3.0 3.5 4.0 4.5 5.0 5.5 6.0 



L [fm] 



L [fm] 



FIG. 6: Relative change in magnetic polarizabilities. Plotted versus L for a few values of the pion 
mass are the finite volume effects A/3a/ for both neutral and charged pions. 



The relative size of finite volume corrections to the magnetic polarizability, ^(3m = 
^"^fa7^" , is investigated in Figured The figure depicts non- negligible finite volume cor- 
rections. The corrections are of comparable size when disconnected diagrams are neglected. 
Carrying out the partially quenched computation, we find 



A r /QTr" 1 connected 

"l[Pm\qem 



25 a 



M JQEM 



f.s. 



4 T 



7/2 



(41) 



IV. SUMMARY AND DISCUSSION 



Above we have investigated finite size effects on pion two-point correlation functions cal- 
culated in constant background fields. Treating the external field strength as weak, couplings 
to pions can be considered as perturbations in coordinate space. One then uses coordinate- 
space perturbation theory on a torus to ascertain the two-point correlation function in finite 
volume. In the case of electric fields, we utilized a time-dependent gauge potential which 
allowed us to sum the couplings to the external field to all orders. The resulting Schwinger 
proper-time Green's functions were used to deduce volume corrections at one-loop order in 
chiral perturbation theory. The effect of such corrections is to renormalize parameters of 
the single-particle effective action in the infrared. To second order in the electric field, we 
found the electric polarizabilities and Born-level charge-squared couplings are renormalized. 
In the case of magnetic fields, we treated the external field perturbatively from the outset, 
and derived the single-particle effective action to second order in the magnetic field. For this 
case, we found the infrared effect on the magnetic polarizabilities, and an identical renor- 
malization of charged-squared couplings as in the electric case. Numerical estimates of the 
finite size effects on polarizabilities show non-negligible dependence on the lattice size for 
pion masses < 350 MeV. Depending on the sign of such corrections relative to the infinite vol- 
ume value, signals for polarizabilities will either be sizably enhanced, or sizably diminished. 
Based on one-loop chiral perturbation theory, the field-squared shifts in the effective mass of 
charged pions will be enhanced in both electric and magnetic fields. On the other hand, the 
shifts for neutral pions will be diminished in both electric and magnetic fields. Results were 
also presented for two-point functions calculated using the lattice approximation of omitting 
the disconnected quark contractions. For the charged and neutral pions this approximation 
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amounts to quenching the electric charges of the sea quarks; while for the neutral pion, the 
annihilation contraction is additionally neglected. 

There are three clarifications that must be made about our results. Firstly and quite 
obviously, the gauge invariance of the single-particle effective actions is not manifest because 
we have expanded results to second order in the external field. In infinite volume, the charge 
couplings are exactly fixed by Ward identities, while in a compact space additional gauge 
invariant terms lead to renormalization of charges and currents. Because we treat the time 
direction as infinite, charge is not renormalized. The possibility to renormalize spatial Born- 
level couplings in a gauge invariant manner has been thoroughly detailed using a momentum- 
space analysis j23|, and we reproduce the findings of that work using a complementary 
coordinate-space method. Finite volume corrections to polarizabilities show up as coefficients 
to field squared terms, £^ and B^. These terms, moreover, can be further split into various 
gauge invariant terms. For example, at second order in the electric field strength one cannot 
distinguish between operators of the form 

Oi = £^ and O2 = 7^Vr(-) ■ £, (42) 

because both terms have the form of polarizability operators at second order. Another way 
to implement a constant electric field on a torus is to use the vector potential 

A^(a;) = (0,0,0,^x3), (43) 

for which the operator O2 vanishes. On a torus, however, the vector potentials, Eqs. ([T]) 
and fH5]) . are not related by a gauge transformation. Consequently the finite volume shift 
of the polarizability (defined as the coefficient of the £^ term of the effective action) will 
be different. A direction for future investigation is to classify all possible gauge invariant, 
single-particle operators allowed on a torus, and determine the running of their coefficients 
with the lattice volume. In this way, one could determine how to implement the external 
field with minimal volume corrections. 

Secondly the finite volume renormalization of charge- squared couplings must be addressed 
to properly fit hadron correlation functions. Let us focus on the electric case. For the charged 
pion, this is relatively straightforward. The two-point function arising from the effective 
action in Eq. fl20l) has precisely the form of Eq. ([5]) projected onto zero three- momentum with 
the mass replaced by the effective mass, namely — >• + A'Ka'^{L) S'^, and the charge- 
squared replaced by the renormalized charge- squared. Because Eq. ([5]) at zero momentum 
is an even function of the charge e, the latter can be achieved simply by the trick e — *■ 
e (1 + |5l[Q^±]). While chiral perturbation theory estimates the latter at the few percent 
level, this expectation can be confirmed by treating the charge-squared as a free parameter 
in fits to charged pion correlation functions.^ The neutral pion, however, presents difficulty 
due to the sign of the renormalized charge-squared coupling. The two-point function derived 



We stress that the charge (calculated from the time component of the current) has not been renormalized, 
as finite volume calculations indeed confirm 23|. For the charged pion with non-zero momentum, the 



spatial current is screened in finite volume, and the effective charge arising from the spatial current is what 
appears in Eq ([5]). The square of this effective charge has been determined above by projecting onto zero 
three-momentum. For the neutral pion, however, charge conjugation invariance forbids renormalization of 
(odd powers of) the current operator (s). In the absence of current screening, however, a renormalization of 
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from Eq. (120!) for the neutral pion appears as Eq. (jSj) projected onto zero three-momentum, 
with — > + 47ra^ (L) and e — > necessitating analytic continuation. The 

two-point function acquires an imaginary part due to the instability in Eq. fl20p for negative 
charge squared. This imaginary part is exponentially small provided either the electric field 
strength is small, E/m^ <^ 1, or the pion Compton wavelength is small compared to the 
lattice size, m^L ^ 1. In practice, this should thus be a very small effect. All neutral 
particles potentially have this difficulty. If the renormalized charge-squared is positive, we 
merely use Eq. ([5]) with Q{L), and their correlation functions should look like the zero- 
momentum projection of charged particle correlation functions. In particular, they will not 
be simple exponentials at long times. If, on the other hand, the renormalized charge-squared 
is negative (as is the case for the neutral pion), our only recourse is to attempt to fit the 
correlation function using the real part of the analytically continued two-point function. This 
procedure relies on the hope that the imaginary part is numerically insignificant to allow 
treatment in Euclidean space. Notice for the connected part of the neutral pion two-point 
function, this issue does not arise because '5l[Q^o]qem^*^'^ = 0- 

Thirdly the choice of coordinate origin is ordinarily arbitrary. In backgournd fields on a 



torus, however, this is no longer the case, as was pointed out for an electric field, see e.g. [19 
Focusing on the electric case and the gauge potential in Eq. ([1]), we observe that two-point 
correlation functions for both charged and neutral particles at finite volume will depend 
on the source time. In infinite volume this was already the case for charged particles via 
Eq. For neutral particles, charge-squared couplings are generated from infrared effects, 
and a neutral hadron source will thus break discrete time-translational invariance. Another 
way to shift the time is through the gauge potential in Eq. ([1]). We are free to translate 
X4 ^ X4 + c in infinite volume as it corresponds to a gauge transformation. In finite volume, 
this shift in time leads to an additional coupling of the matter fields to a constant gauge 
potential 

Al^\x) = {0,0,-Sc,0), (44) 

that cannot be gauged away on a torus. The effect is equivalent to a twisted boundary 
condition on the matter field, see e.g. [411]. Consequently a hadron with charge Q acquires 
kinematic momentum, k = —QScz. One can extend the analysis formulated here to this 
flavor-twisted case by working at non-zero spatial momentum. Volume effects can be com- 
puted [i^, however, the expectation is that broken symmetries will lead to extra finite 



volume terms [43|, 14^, |45|. This complication notwithstanding, effects on the two-point 
function can be determined. On general grounds, we expect all terms in the single particle 
effective action to be renormalized including the current, or equivalently the momentum. 
As was the case for periodic boundary conditions, neutral particle correlation functions will 
no longer have a simple exponential falloff at long times. 

Using the methods established here, many subsequent investigations are possible. Of 
particular interest is the calculation of finite volume effects for baryon electromagnetic prop- 



the Thomson cross section gives rise to a charge-squared coupUng. This is permitted because two current 
operators are involved. The finite volume neutral pion propagator at non-zero momentum is then given 

by 

where (2:4, s| 2:4, 0) depends on the effective charge squared (5l[Q^o]- 
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erties. These too can be addressed using the present framework. In order to rehably extract 
these quantities from lattice QCD simulations in background fields, one must ascertain 
whether volume corrections are as sizable as in the meson sector. 
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